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On 

Abstract 

The fluctuations of the vacuum energy are treated as a non-equilibrium pro- 
cess and a stochastic model for the cosmological constant is presented, which 
yields a natural explanation for the smallness or zero value of the constant 
■ in the present epoch and its large value in an era of inflation in the early 
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0^ ■ universe. 
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It is widely considered that the problem of the cosmo logical constant has not been 
satisfactorily explained. The problem is this: the vacuum fluctuations of particle fields add 
up to a contribution which is much larger than the observational limits by many orders of 
magnitude. Attempts using quantum cosmology, based on e.g., Euclidean path integrals and 
wormholes, have been criticized for various technical reasons, which question it as a possible 
solution to the problem [|J. Unbroken supersymmetry can "protect" the constant A from 
becoming non-zero, but supersymmetry is badly broken in the real universe at an energy 
> 1 TeV. No other symmetry is known to exist in Nature that guarantees that A remains 
zero. Recently, a solution using the Jordan, Brans-Dicke theory of gravitation was proposed 
in which the cosmological "constant" is a smooth function of time 0. 

It is believed by many that the problem would be solved within a physically consistent 
theory of quantum gravity. Since such a theory does not presently exist, we shall propose a 
possible solution to the problem based on treating the vacuum energy as a stochastic system. 
We describe this system by a phenomenological equation for the cosmological constant A, 
which appears in Einstein's field equations: 

R^u - -xg^R + Asv = (1) 
The cosmological constant enters through the vacuum energy density: 

Tvtiv = —PvQiiv = ~^q3^v (2) 

Today, A has the incredibly small value, A < 10 -46 GeV 4 , whereas generic inflation models 
require that A has a large value during the inflationary epoch. This is the source of the 
cosmological constant problem. 

Let us treat the vacuum energy as a fluctuating environment and consider A as a variable 
characterizing the state of this system, which is determined by the equation : 

A = aA - A 2 , (3) 

where A = dA/dt, a is a parameter which corresponds to the difference between the rate of 
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creation and annihilation of particles in the vacuum. The second term is a self-restriction 
term which limits the growth of A. The solution of ([3]) is given by 

A(t) = A(0)exp(at){l + A(0)[(exp(at) - l)/ar]} -1 . (4) 

At a = 0, there is only one stable stationary state solution A = 0, and at this point the 
solution becomes unstable and bifurcates into a new branch of steady stable solutions, A = a. 
At a = the system undergoes a second-order phase transition. 

Let us consider the situation in which the vacuum fluctuations are rapid compared with 
Tmacro = ct' 1 , which defines the macroscopic scale of time evolution. We shall assume that 
the parameter a can be written as a t = a + cr£ t , in which a is the average value, £ t is 
Gaussian noise and a measures the intensity of the vacuum fluctuations. Then, we write 
Eq.(|) as 

dA t = (aA t - A 2 t )dt + aA t dW t = f(A t ) + ag(A t )dW t , (5) 

where dW t is a Wiener process. We shall use the Ito integral to describe the diffusion 
process, although the Stratonovich integral would predict the same qualitative results ||. 
The probability density p(x, t) satisfies the Fokker-Planck equation: 

a 2 

d t p(x,t) = -d x [(ax - x 2 )p(x,t)} + —d xx {x 2 p(x,t)). (6) 

The diffusion process is restricted to the positive real half line and and oo are intrinsic 
boundaries, because <?(0) = and /(oo) = — oo. The probability of the diffusion process 
reaching infinity as t — > oo is zero, since infinity is a natural boundary. Moreover, zero 
is a natural boundary if a > a 2 /2 0, so neither boundary is accessible and no boundary 
conditions need be imposed on the Fokker-Planck equation. For a < cr 2 /2, it can be shown 
that zero is an attracting boundary. 

The stationary-state solution for the probability density, p s (x), of Eq.(|6|) is given by 

p s {x) = Nx( 2a ^- 2 e W (-^p). (7) 
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The normalization constant N is 




(8) 



If p(x,t) is integrable between and oo, then a stationary state solution exists when a > 
cr 2 /2. If it does not exist, then the probability density will be concentrated at zero, i.e., 



The extrema of p s play the role of order parameters for non-equilibrium phase transitions, 
and they are determined by the equation: 



They are given by A ml = and A m2 = a — cr 2 (a > a 2 ). The maximum A m2 always exists, 
while the maximum A ml exists for < a < cr 2 . Two transition points exist: one at a = a 2 /2 
and one at a = a 2 . We therefore have the following situation: (a) If a < a 2 / 2, then zero is 
a stable stationary point for A. (b) The point a = a 2 /2 is a transition point, since A = 
becomes unstable and a new stationary probability density is produced, (c) The stationary 
density becomes divergent at A = when a 2 / 2 < a < a 2 . Although zero is no longer a 
stable stationary point, it remains the most probable value, (d) For a = a 2 , the probability 
density p s {A) undergoes a transition, in which A can take on large values by increasing the 
intensity of the fluctuations, even though the average state of the vacuum is kept constant. 

For < a < cr 2 /2, the vacuum fluctuations dominate over the growth or decline of A, 
although the value zero is still the most probable value for A, since the distribution function 
has a vertical slope at A = 0. Because we are using a continuous variable, A never reaches 
the boundary zero in a finite time. This model breaks down when A is vanishingly small 
and the probability of having A = is defined for < e << 1 and 0<a:<cr 2 /2by[|5|: 



When a > a 2 /2, the growth of A dominates the influence of the vacuum fluctuations, 
and in the neighborhood of zero the probability of A = drops to zero. At a = a 2 /2 real 
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p s (A) = 5(A) for a < a 2 /2. 



A 2 m -(a- a 2 )A m = 0. 



(9) 




(10) 



growth of A becomes possible corresponding to a change from a degenerate random variable 
for steady-state behavior to a stochastic variable; the boundary at a = switches from 
attracting to natural. The transition point a = a 2 corresponds to a qualitative change in 
the stochastic variable A with no change in the nature of the boundary. The probability of 
A = drops abruptly to zero. This phenomenon is inherently nonlinear. 

We now see the following scenario emerging from our model. In the inflation era, the 
intensity of vacuum fluctuations is large and a > a 2 , causing a second-order phase transition 
and a maximum in A not near zero. This corresponds to the large vacuum energy needed to 
drive inflation [j/j. As the universe expands the intensity of vacuum fluctuations decreases 
and for < a < a 2 /2 or a 2 /2 < a < o 2 the probability density is largest when A is non- 
vanishing and small, which can lead to a current value of A that can be used to fit the 
observational data. If the stationary probability density p s does not exist for a < cr 2 /2, then 
A = is a stationary point; the drift and diffusion vanish simultaneously for A = and 
p s (A) = 5(A). This corresponds to the case when A is vanishingly small. 

Thus, our model provides a natural explanation, in terms of non-equilibrium stochastic 
processes in an expanding universe, for the behavior of A needed to fit observational data 
and still be consistent with inflationary models. 

According to general relativity, the equation that governs the expansion factor R(t) is 
given by 

where pu is the mass density; k = —1, 0, +1 and H is the Hubble constant, whose observable 
value at present time to is denoted by H . We define 

n tot = n M + n A = i-n k . (12) 

Many observational tests can constrain but the gravitational lensing method provides 
the most direct constraint with the result: Qq,a < 0.75 |J. A recent analysis of the cosmo- 
logical data showed that for f^A = 0.65 ± 0.1, VLm — 1 — and a small tilt: 0.8 < n < 1.2, 
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models exist which are consistent with the available data and an inflationary spatially flat 
universe 0. 
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